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Rare-earth-based permanent-magnet materials rich in iron have relatively low ferromagnetic or- 
dering temperatures. This is believed to be due to the presence of antiferromagnetic exchange 
interactions, besides the ferromagnetic interactions responsible for the magnetic order. The mag- 
netic properties of Ce2Fei7 are anomalous. Instead of ferromagnetic, it is antiferromagnetic, and 
instead of one ordering temperature, it shows two, at the Neel temperature Tn ~ 208 K and at 
T t ~ 124 K. Ce2Fei7, doped by 0.5% Ta, also shows two ordering temperatures, one to an antifer- 
romagnetic phase, at Tn ~ 214 K, and one to a ferromagnetic phase, at To ~ 75 K. In order to 
clarify this behavior, single-crystalline samples were prepared by solution growth, and characterized 
by electron microscopy, single crystal x-ray diffraction, temperature-dependent specific heat, and 
magnetic field and temperature-dependent electrical resistivity and magnetization. From these mea- 
surements, magnetic H-T phase diagrams were determined for both Ta-doped Ce2Fei7 and undoped 
Ce2Fei7. These phase diagrams can be very well described in terms of a theory that gives magnetic 
phase diagrams of systems with competing antiferro- and ferromagnetism. 

PACS numbers: 75.30.Kz,64.60.Kw,71.20.Lp,61.50.Nw 
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I. INTRODUCTION 

Modern permanent magnet materials are intermetal- 
lic compounds containing rare earths (R) and tran- 
sition metals, which order magnetically at temper- 
atures much above room temperature, and combine 
a large ferromagnetic moment with a large easy-axis 
magnetic anisotropy 1 . The large magnetic moment 
(and the high ordering temperature) is mainly provided 
by strongly (and ferromagnetically) coupled itinerant 
transition-metal magnetic moments. The strong mag- 
netic anisotropy is mainly provided by the localized 4f- 
rare-earth magnetic moments, to which the transition- 
metal magnetic moments are coupled. 

Because of its relative abundance in Earth's crust, Fe is 
the most favorable (cheapest) magnetic transition metal 
for use in these permanent-magnet materials. Among the 
rare-earth elements, Ce is most abundant, and therefore 
cheapest, potentially making an Fe-rich Ce-Fe compound 
an economical permanent magnet material. However, the 
aforementioned requirements are not met. Instead, the 
Ce-Fe compound richest in Fe, rhombohedral Tli2Zni7- 
type Ce2Fei7 shows abnormal- magnetic behavior. How- 
ever, since Ce2Fei 7 is chemically closely related to suc- 
cessfully applied permanent magnet materials^, such as 
SmCo5 and Nd2Fei4B, understanding its magnetism may 
lead to a greater understanding of magnetism in rare- 
earth-transition-metal intermetallic compounds. 



The R2Fei7 compounds form, dependent on the size 
of the rare-earth ion, in two different but related crys- 
tallographic structures, rhombohedral Tl^Zniytype for 
light (larger) rare earths and hexagonal Tli2Nii7-type for 
heavy (smaller) rare earths 4 . Generally, they are ferro- 
magnetic (Fe moments parallel to R-4f moments) for light 
rare earths, and ferrimagnetic (Fe moments antiparallel 
to R-4f moments) for heavy rare earths^. Exceptions are 
Lu2Fei7 and Ce2Fei7. These compounds have been re- 
ported to show, at least in some temperature range, an- 
tiferromagnetic behavior (see e. g. Refs. 0(6.7), which for 
Lu2Fei7 occurs between the magnetic ordering tempera- 
ture (275 K) and a second critical temperature (140 K), 
below which it is ferromagnetic. Under a hydrostatic 
pressure of 0.4 GPa, Lu2Fei7 remains antiferromagnetic 
at temperatures down to 5 K£. 

As a representative for the magnetic behavior of the 
Fe-magnetic subsystem in R-Fe compounds, one often 
considers the Y-Fe analogs (La-Fe binary intermetallic 
compounds do not form), and such an analog for Ce2Fei 7 
is Y 2 Fei 7 . Y 2 Fei7 can be prepared either with a hexago- 
nal Tli2Nii7-type or with a rhombohedral Tli2Zni7-type 
crystal structure. Hexagonal and rhombohedral Y2Fei 7 
are very similar in their magnetic properties, but hexag- 
onal Y2Fei 7 may be easier to obtain (see e. g. Ref. M), 
while rhombohedral Y2Fei 7 is easier to modeU^, due to 
its simpler crystal structure. Hexagonal Th 2 Nii7-type 
Y2Fei7 is an easy-plane ferromagnet: it orders magneti- 
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cally near 310 K£>£, and the saturated magnetic moment 
equals ~ 1.95 ^e/Fe at 4.2 K. In contrast to the mag- 
netic behavior of Co-rich Y-Co compounds, Fe-rich Y-Fe 
compounds have magnetic ordering temperatures that 
decrease with increasing Fe concentration, which helps 
put the relatively low ordering temperature of Y2Fei7 
in perspective (see e. g. Ref. y). Electronic-structure 
calculations indicate that Y-Fe compounds can be clas- 
sified as weak ferromagnetsi -, with incomplete occupa- 
tion of the majority-spin bands. For Y 2 Fe!7, a rapid 
variation of the density of states near the Fermi levels 
was found. This may explain the calculated strong (and 
non-linear) volume dependence of magnetization. Exper- 
imentally, hexagonal Tli2Nii7-type Y2Fei7 was found to 
have a negative thermal expansion below its magnetic 
ordering temperature^, which was phenomenologically 
related to strongly distance-dependent positive and nega- 
tive magnetic exchange interactions between different Fe- 
magnetic moments on different crystallographic sites, a 
notion that was confirmed by calculations related to elec- 
tronic structure^ 2 -. Recently, Prokhnenko et al. — found 
that hexagonal Y2Fei7 shows antiferromagnetic behavior 
under sufficiently high pressures. They indicate that the 
magnetism in Y2Fei7 under pressure may be compared 
to antiferromagnetism in unstable 7-Fe, which is still a 
challenge for electronic structure calculation s 13 ' 14 ' 15 ' 16 . 

As determined from polycrystalline samples Ce2Fei7 
orders antiferromagnetically below Xn ~ 208 K — . At 
125 K, a second magnetic-order transition is observed, to 
another cmiiferromagnetic phase £ i 17 i 18 i 19 i 20 i 21 i 22 . How- 
ever, there arc reports, in which Ce2Fei7 is found to 
be /erromagnetic at temperatures below an antiferro- 
ferromagnetic transition temperatur o 21 ' 22 ' 23 ' 24 , which is 
sample dependent and ranges between ~ 20 and ~ 
120 K 2 — . It has also been reported^ 5 -, that Ce2Fei7 doped 
with ~1.4% Si is ferromagnetic below ~ 120 K, and an- 
tiferromagnetic between ~ 120 K and ~ 210 K. Based 
on their Mossbauer spectroscopy experiments on differ- 
ent samples, nominally undoped or doped with Si or Al, 
from different laboratories, Hautot et al. — found that 
the samples which are ferromagnetic at low temperatures 
may be doped by small amounts of Si or by Al. Such 
doping may be caused by samples reacting with AI2O3 
crucibles 3 -^ or silica ampoules. 

Based on the lattice constants, the unit cell of 
Ce2Fei7 is smaller than expected from La-contraction, 
Buschow and van Wieringen 2 - deduced that Ce in Ce2Fei7 
is tetravalent, and therefore carries no local 4f mo- 
ment. Later X-ray absorption spectroscopy, by Isnard et 
a ^ 28, 29,30 an( j Vandormael et al^ indicated that Ce in 
Ce2Fei7 is in a mixed-valent state. Furthermore in sub- 
stitutionall y 30 ' 31 and interstitiall y 28 i 29 modified Ce2Fei7, 
the spectroscopic Ce valence decreases with increasing Ce 
site volume. 

Powder-neutron diffraction experiments on apparently 
undoped Ce 2 Fei 7 i£i2 2 - indicated that below T N - 208 K, 
Ce2Fei7 orders in a complex AF structure, which is mod- 
ified below T t ~ 125 K. In both cases the magnetic mo- 



ments were found to order perpendicular to the (hexag- 
onal notation) c-axis. These powder-neutron results will 
be discussed elsewhere in this volume 3 -^. 

Magnetization experiments on polycrystalline samples 
hinted at a rich magnetic (H-T) phase diagram 6 -^, that 
is strongly pressure dependent—, and involves strong 
magnetovolume effects 21 . Recently, Makihara et al. — 
succeeded in preparing single crystals of apparently un- 
doped Ce2Fei7, that displayed anisotropic magnetiza- 
tion. However, a magnetic phase diagram of single- 
crystalline Ce2Fei7 has not been determined. 

Here we report on solution-grown single crystals of 
Ce2Fei7 and of Ta-doped Ce2Fei7, and on the character- 
ization of the crystals by magnetization, electrical trans- 
port, specific heat, x-ray diffraction, and electron-probe 
micro analysis. Below, we will first study the composi- 
tion and the crystal structure of Ta-doped crystals and 
compare them to undoped crystals. Then we present the 
results of extensive thermodynamic and transport mea- 
surements to determine the anisotropic H-T phase dia- 
grams for undoped Ce2Fei7. After this, we present an 
H-T phase diagram obtained on a Ta-doped crystal of 
Ce2Fei7. Finally, the H-T phase diagrams are discussed 
in terms of a theory of competing ferro- and antiferro- 
magnetic interactions. 



II. EXPERIMENTAL 

Traditionally, R-Fe compounds are grown by self-flux 
out of Ta crucible s 33 ' 34 . This technique has been used 
to grow Nd2Fei7 as well as Nd2Fei4B single crystal s 35 ' 36 . 
Unfortunately, as will be discussed below, Ce2Fei7 allows 
for a slight Ta uptake^ 2 - on one of the Fe-crystallographic 
sites, which profoundly alters the magnetic behavior of 
the sample. In order to fully control the amount of Ta in 
the sample, we have used MgO crucibles (Ozark Techni- 
cal Ceramics, Inc.), as described in Ref.[22l to grow single 
crystals of Ce2Fei7 as well as of intentionally Ta-doped 
Ce 2 Fei 7 . 

Starting alloys were made by arc melting pieces of 
Fe(99.98%, SCM Metal Products Inc., Cleveland, OH), 
Ce (99.9% elemental, 99.98% metals basis, Materials 
Preparation Center, Ames Laboratory—) and Ta (99.9%) 
in a standard arc furnace under a partial pressure of ~ 0.8 
bar of Ti-gettered ultra-high-purity Ar. Prior to the crys- 
tal growth experiments, the composition and tempera- 
ture range for crystal growth were optimized with the 
aid of differential thermal analysis 2 ^ (DTA) . 

The button of the starting alloy was placed in the 
growth crucible, a 2 ml MgO crucible. This growth cru- 
cible was capped with a strainer, a piece of Ta foil with 
several holes in it, while care was taken that the foil did 
not come in contact with the alloy, to prevent possible 
dissolution of the Ta. On top of the capped MgO cru- 
cible, a catch crucible, an inverted 2 ml AI2O3 crucible, 
was placed. This assembly was sealed in an amorphous 
silica ampoule under ~ 0.3 bar of ultra-high-purity Ar- 
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gon. 

To obtain undoped Ce2Fei7, an alloy with composi- 
tion Ceo.475Feo.525 was heated to 1200°C, then allowed to 
equilibrate for ~ 2h, quickly cooled to 1000°C, and then 
cooled slowly to 940°C, where the flux was decanted. The 
liquid alloy had visibly reacted with the crucible: the cru- 
cible had become brown, and on the bottom, apparently 
up to the level of the liquid, a brownish crust had formed. 
On this crust, which contained Ce oxide, we also found 
crystals of Ce2Fei7, which we removed mechanically. As 
we will see below, no Mg was found in these crystals, as 
may be expected from the limited solid solubility of Mg 
in Fe 3 ^. 

Crystals typically grow blocky-rhombohedral or tab- 
ularly, with clean and clear facets, see Fig. [TJ For the 
tabular crystals, the facets with largest surface area were 
found perpendicular to the hexagonal c-axis, as con- 
firmed by x-ray diffraction. (From here, we will use 
hexagonal notation for the rhombohedral structure.) It 
appeared that the probability to obtain tabular (or even 
plate-like) crystals is larger when the alloy was cooled 
more quickly, and conversely, the probability to obtain 
blocky (more 3-D) crystals is larger when the alloy was 
cooled more slowly. Therefore, cooling rates varying be- 
tween 0.3 - 2°C/h were used. (Note that the DTA- 
optimization^ 7 - is useful here: slow cooling was performed 
over only 60° C.) Among the blocky crystals we could not 
find crystals longer, in the c-direction, than ~ 1.5 mm. 
Moreover, crystals reasonably long in the c-direction ap- 
peared to have inclusions that became visible when the 
crystals were cut. 

Quantitative electron-probe microanalysis (EPMA) 
was performed with a JEOL JXA-8200 Superprobe with 
a 20 kV acceleration potential and a 20 nA beam current. 
Crystals of undoped Ce2Fei7 and of Ta-doped Ce2Fei7 
were polished and mounted side by side. As a stan- 
dard for Ce, a single crystal of the air-insensitive binary 
line compound CeRu2 was used. For Mg, a sample of 
Mg2Si04 was used as a standard, and for Fe and Ta, 
pieces of those elements were used as standards. 

Single-crystal x-ray diffraction measurements were per- 
formed on the undoped and Ta-doped Ce2Fei7. Crystals 
with diameter ~ 10 /zm 3 were obtained by crushing crys- 
tals from the respective batches. Room-temperature x- 
ray diffraction data were collected on a STOE IPDSII 
image plate diffractometer with Mo Ka radiation, and 
were recorded by taking 1° scans in 9 in the full recipro- 
cal sphere. 29 ranged from 6° to 63°. Numerical absorp- 
tion corrections for both crystals were based on crystal 
face indexing, followed by a crystal-shape optimization. 
Structure solution and refinement were done using the 
SHELXTL program. 

Magnetization measurements at temperatures between 
2 K and 300 K in applied fields of up to 5.5 or 7 T were 
performed in Quantum Design MPMS-5 and MPMS-7 
magnetometers, respectively, on plate- like crystals as well 
as on blocky crystals. Zero-field ac-susceptibility was 
measured in a Quantum Design MPMS-5 magnetome- 



ter, excitation frequency 1000 Hz, amplitude 0.2 mT, on 
a tabular crystal of ~ 97 mg. Specific heat was deter- 
mined in a Quantum Design PPMS system on the same 
crystal as the one used for ac-susceptibility. For electrical 
transport measurements, we used a standard four-probe 
technique. Bars were cut with a wire saw, both for the 
current j flowing parallel to the c-axis and for the cur- 
rent j flowing perpendicular to the c-axis, and electrical 
contact was made with Epo-tek H20E silver epoxy, with 
typical contact resistances of about 4 f2. To obtain a bar 
for j_L c-axis was relatively easy: a plate-like crystal ~ 
0.15 mm thick, and several mm diameter was used to cut 
a bar of 3 mm length and 1 mm width. To obtain a good 
sample for j / / c-axis was much more difficult, because 
of the limited available length. The resulting sample was 
0.4 mm thick, 0.7 mm wide and 1.1 mm long. 



III. CHEMICAL AND STRUCTURAL 
CHARACTERIZATION 

A total of four single-crystal samples were examined 
by EPMA. For one of these, the starting alloy nominally 
consisted only of Ce and Fe, and for the other three, 
the starting alloy was doped with Ta (by 0.05%, 0.5% 
and 1%). The nominally undoped sample was very care- 
fully checked for both Mg and for Ta by performing slow 
wavelength-dispersive spectroscopy scans over the Mg-K 
and Ta-L lines, respectively, on the standards as well as 
on the sample. The spectra from the nominally undoped 
sample showed no discernible Mg or Ta peaks. Spectra 
from Ta-doped samples did produce Ta peaks, and Ta- 
peak heights obtained by scanning were consistent with 
Ta-peak heights obtained by a three-point analysis. The 
analysis indicated that the undoped sample contained 
less than ~ 300 ppm Mg and less than ~ 300 ppm Ta. 
Quantitatively, the composition of the undoped sample 
was found equal to C^Feie.gfi)- The nominally 0.05%, 
0.5%, and 1% Ta-doped samples were found to contain 
0.42(5)%, 0.37(5)%, and 0.57(4)% Ta, respectively, which 
indicates that there is a solid-solubility maximum for Ta 
in Ce2Fei7. Quantitatively, the compositions were found 
equal to Ce2Fe 17 . 1( 2)Ta .o 8 (i), Ce2Fe 17 .o(i)Ta .o7(i), and 
Ce2Fe 17 . oo(9)Ta . 108(7) j f° r nominally 0.05%, 0.5% and 
1% Ta, respectively. Note that these compositions do 
not differ significantly, i.e. by more than two standard 
deviations, from (2-17) stoichiometry. EPMA also in- 
dicated that samples may have parasitic inclusions of 
CeFe2, which was also noted in Ref . l2~il. 

An undoped crystal and the EPMA-crystal, grown out 
of the nominally 0.05% Ta-doped alloy were used for 
single-crystal X-ray diffraction experiments. The crys- 
tallographic and structural data are summarized in Ta- 
bles U and [Til for undoped and Ta-doped Ce2Fei 7 , re- 
spectively. We refined the x-ray diffraction-data for the 
undoped sample assuming only Ce and Fe form the struc- 
ture. The diffraction data indicates the crystal has the 
rhombohedral Tli2Zni7-type crystal structure with space 
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Atom(site) 


x/a 


y/b 


z/c occ. U eq 


Ce(6c) 








0.3436(1) 1 0.0068(3) 


Fe(6c) 








0.0965(2) 0.954(9) 0.0074(11) 


Ta(6c) 








0.0965(2) 0.046(9) 0.0074(11) 


Fe(9d) 


1 

2 





| 1 0.0060(5) 


Fe(18/) 


0.2915(2) 





1 0.0071(4) 


Fe(18/i) 


0.5012(1) 0.4988(1) 0.1549(1) 1 0.0080(4) 





E 



Ce(6c> 

^ Fa (6tt 

J Fa (9d) 

4 Fe (16f) 

• Fb(1«i) 



FIG. 1: Schematic drawing*^ of the rhombohedral unit cell 
of Ce2Fei7, with hexagonal axes. In Ta-doped Ce2Fei7 the 
Fe (6c) site is partially occupied by Ta. The inset shows a 
photograph of a plate-like crystal on a mm-grid background. 
The c-direction is perpendicular to this grid, and the sides 
coincide with main planar directions. 
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0.0086(2) 


Fe(18/i) 


0.50150(4) 0.49850(4) 0.15502(5) 
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TABLE I: Atomic coordinates and equivalent isotropic dis- 
placement parameters (A 2 ) for Ce2Fei7. \J eq is defined as 
one third of the trace of the orthogonal U !J tensor. Space 
group R3m, a = 8.4890(12) A, c = 12.410(4) A, R = 0.0149, 
R w = 0.0361. 



group i?3m, as shown in Fig. [T] with a full occupancy of 
all available sites. The high quality of the refinement (the 
refinement factor R = 0.0149), indicates that this crys- 
tal has no obverse-reverse twins (reticular merohedry) , as 
found for the related compound Pr 2 Fei7, see Ref. Ho. We 
found for the lattice parameters of the undoped Ce2Fei7 
crystal a = 8.4890(12) A, and c = 12.410(4) A. These 
unit-cell parameters agree very well with those found for 
apparently undoped Ce2Fei7, see Ref. [U The atomic 
position parameters, see Table UJ are in good agreement 



TABLE II: Atomic coordinates and equivalent isotropic dis- 
placement parameters (A 2 ) for Ta-doped Ce2Fei7. XJ eq is de- 
fined as one third of the trace of the orthogonal U !J tensor. 
Space group R3m, a = 8.4972(12) A, c = 12.447(3) A, R = 
0.0496, R w = 0.0535. 



with published values^. 

Although it was clear from the EPMA experiments 
that there was Ta in the doped sample, we initially re- 
fined the diffraction data from this crystal assuming that 
only Ce and Fe form the structure. Although in this case 
the i?-factor was only slightly higher than the ii-factor 
for the final refinement shown in Table ITU 0.0504 instead 
of 0.0496, the refined isotropic displacement parameter 
for the Ye-6c site became very small, 0.0021(6) instead 
of 0.0074(11). This value is much lower than the dis- 
placement parameters found for the other Fe atoms in 
the structure, whereas for pure Ce2Fei7, Table U values 
were about equal for the different Fe sites. Since small 
thermal displacements have the same effect on diffracted 
beam intensities as larger electron densities, it appears 
that the electron density on the Ye- 6c site is too low in 
the model, without Ta. The introduction of an Fe/Ta 
mixture on the Ye- 6c site brings the temperature factor 
to the level of the other Fe atoms. Further refinements, 
allowing for mixtures on the other crystallographic sites, 
made clear that Ta has a strong preference to substitute 
for Fe on this Ye- 6c site. Therefore, the structure was 
refined with Ta only on that site. The final refinement 
factor R — 0.0496, indicates a good fit. Residual electron 
density not described by the current structural model 
does not lead to Ta occupying possible interstitial sites by 
more than about 0.5%. The refined amount of Ta leads 
to an empirical x-ray formula of Ce2Fei6.gi(2)Tao.og(2) , 
which is, within experimental error, the same as the com- 
position found by the EPMA experiment for this crystal, 

Ce2Fei 7 .i(2)Ta o .08(i)- 

The lattice parameters and unit-cell parameters for 
the Ta-doped Ce2Fei7 crystal: space group i?3m, a = 
8.4972(12) A, c = 12.447(3) A, are slightly larger than the 
lattice parameters for the undoped crystal, which is con- 
sistent with reported 21 lattice parameters for Ta-crucible 
grown Ce2Fei7. 

The result that Ta substitutes for Fe on the Ye- 6c site, 
is consistent with expectations, and the results that Ta 
has a solid-solubility maximum in Ce2Fei7, is both con- 
sistent with expectations and plausible. All transition 
metals with a radius larger than Fe, which includes Ta, 
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Ref. H3, prefer to substitute for Fe on the Fe-6c site^ 3 - 
in Nd 2 Fei7 with the same crystal structure as Ce 2 Fei7. 
The solid solubility of Ta in isostructural Sm 2 Fei7 has^ 
a maximum of 2.3 at% of Ta in Sm2Fen- x Ta x . In solu- 
tion growt h 27 ' 33 i 34 , the composition of the product (crys- 
tals) is different than the composition of the initial melt 
(crystals plus solvent). Here, the mass ratio (yield) of 
crystals-to-flux was typically about 10%, which makes it 
theoretically possible to obtain crystals with as high as 
0.5 at% Ta for an initial melt with 0.05 at% Ta. 

The three Ta-doped crystals under investigation are all 
doped by about 0.5 % Ta, which indicates that there ex- 
ists, besides a limited solid solubility of Ta in Ce2Fei7, a 
mechanism that favors Ta to dope the crystals (Ce2Fei7) 
rather than remain in the solvent (~ Ceo.sFeo.s). This 
mechanism might be related to the limited solubility of 
Ta in liquid Ce (~ 0.005 % at - 1000°C), see Ref. 0. 



IV. PURE CE 2 FEi 7 - ZERO-FIELD RESULTS 

Temperature-dependent specific heat of Ce2Fei7 is pre- 
sented in Fig. Two A-type peaks, at 124 K and at 
209 K, are clearly seen. These peaks are consistent with 
magnetic ordering transitions found earlier in polycrys- 
talline samples^. 

Results of anisotropic zero-field ac-susceptibility mea- 
surements for both H ac _L c-axis and H ac / / c-axis, are 
presented in Fig. ^p. For both ac-field directions, three 
features characterize the data: a sharp rise near 225 K, 
and two weaker anomalies, near 208 K and 125 K, re- 
spectively. The temperature of the sharp rise, ~ 225 K, 
is consistent with the ferromagnetic ordering tempera- 
ture of ~ 230 K of CeFez^- (Inclusions of CeFe2 were ob- 
served in samples used for EPMA, see above.) Moreover, 
the sample used for the ac-susceptibility presented here 
was the same as the one used for the specific heat exper- 
iments described above. In that experiment, no anomaly 
was observed near 225 K. From a magnetization isotherm 
measured at 5 K (not shown) we estimate the amount of 
CeFe 2 in this sample as less than 0.5%. As for the other 
two anomalies, they occur at temperatures near the tem- 
peratures of the peaks in specific heat, and are consistent 
with the previously reported ordering temperatures for 
Ce 2 Fei 7 . 

Results of zero-field resistivity measurements are pre- 
sented in Fig. [2]2. Because of the small dimensions of 
the samples, the estimated error in the determination of 
the resistivity was relatively large, about 20 %. At room 
temperature, above the magnetic phase transition tem- 
peratures, overlapping values were found, for j / j c-axis 
and for j _!_ c-axis, 150 (30) /xf2 cm and 190 (38)/if2 cm, 
respectively. Because of this overlap, the existence of 
anisotropy in the resistivity at room temperature is not 
clear. To further study possible anisotropic temperature 
dependence, the resistivity results were normalized to the 
300 K values. 

In polycrystalline samples, resistivity anomalies asso- 
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FIG. 2: a Zero-field specific heat C p as a function of temper- 
ature. The A-type peaks, at T t = 124(1) and T N = 209(2) K, 
are due phase transitions, as described in the text, b ac- 
susceptibility for the ac-field applied both _L c-axis and // 
c-axis. The sharp rise, in both curves, near 230 K is due to 
parasitic CeFe2. The dashed lines indicate the temperatures 
where peaks were found in the specific heat, c Normalized 
temperature-dependent zero-field resistivity p(T) for the cur- 
rent flowing parallel to the c-axis (open circles), and for the 
current flowing perpendicular to the c-axis (closed circles). 
The dashed lines indicate the temperatures where peaks were 
found in the specific heat. 



ciated with the magnetic ordering near 209 K were not 
observe d 6 i 18 ' 31 . Here, above 220 K, the behavior of the 
resistivity appears indistinguishable for both employed 
current directions. Below that temperature, thus near 
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209 K, the temperature of the highest-temperature tran- 
sition found in specific heat, the normalized zero-field 
resistivity for j _L c-axis is lower than for j // c-axis. 
Thus near the magnetic ordering temperature, the resis- 
tivity behavior becomes anisotropic. Both curves appear 
approximately linear with temperature between ~ 130 K 
and - 190 K. 

At a still lower temperature of ~ 124 K, consistent with 
the A- type peak in specific heat, an anomalous increase 
of the resistivity occurs, for both current directions, sim- 
ilar to the resistivity behavior below the ordering tem- 
peratures of heavy rare earths like Er or Tm (see e. g. 
Refs. Ho] and [I?], respectively). Such an increase, though 
much less pronounced, was observed earlier in polycrys- 
talline Ce2Fei7^^. The increase is much larger for j / / 
c-axis than for j _!_ c-axis. This is consistent with the 
formation of anisotropic superzone gaps, gapping parts 
of the Fermi surface, particularly _L c-axis. 

We estimate the gapped fraction of the Fermi surface 
by considering the conductivity 4 ^ 



12n 3 h 



AdSv 



(1) 



with A the electron mean free path, and the integral over 
the Fermi surface area Sf- Assuming that A does not 
change anomalously due to the transition, the gapped 
fraction of the Fermi surface is proportional to the rela- 
tive reduction 6 in conductivity a. If we further assume 
that all parts of the Fermi surface contribute equally to 
conductivity, the gapped fraction is equal to the relative 
reduction in conductivity. 

Then the assumed resistivity without superzone gaps, 
Pwithout, just below 125 K , follows the apparently linear 
behavior observed above 125 K. This resistivity without 
superzone gaps below 125 K is obtained by linear extrap- 
olation of the resistivity between ~ 130 and 180 K. An 
estimate for 5 is obtained from Eq. [1] 



6 = 1- 



^without 
Pwith 



(2) 



The reduced resistivity for j / / c-axis shows a maxi- 
mum near 107 K, and we assume for our estimate that 
the gaps have fully developed there. The reduced resis- 
tivity at 107 K reads p w ith,rcd=0.92, and the extrapolated 
reduced resistivity reads /O w ;thout,red=0.83. Therefore, we 
estimate that the conductivity for j / / c-axis is reduced 
by 5 ~ 10 % due to the superzone gaps. Obtained in a 
similar manner, we estimate that the conductivity for j 
_L c-axis is reduced by ~ 3%. This anisotropic behav- 
ior is consistent with expectations due to the reported 
propagation vectors for the magnetic structure, parallel 
to the c-axi a 1:L i 22 i 49 . 

Furthermore, the zero-field resistivity for both current 
directions shows a rather small residual resistance ratio 
(RRR), about 2 for j j j c-axis and about 3 for j _L c- 
axis. Such values may indicate poor sample quality, how- 
ever, because in an applied field RRR is much larger for 



both current directions, as will be shown below, we think 
that in Ce2Fei7 a strong magnetic scattering is the major 
cause for such low RRR values. 



■2 : j // c-axis 




j // c-axis 



T(K) 
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FIG. 3: Left: a Specific heat, b 8{xT)/dT for H ac ± c-axis 
(closed circles) and for H ac / / c-axis (open circles) , and c 
dp/dT for the current j / / c-axis (open circles), and for j 
_L c-axis (closed circles), between 110 and 135 K. Right: d 
Specific heat, e d(xT)/8T for H ac _L c-axis (closed circles) and 
for H ac // c-axis (open circles) , and f dp/dT for the current 
j /j c-axis (open circles), and for j _L c-axis (closed circles), 
between 195 and 220 K. 



As already indicated above, in the measurements of 
specific heat, ac-susceptibility and resistivity, the two or- 
dering temperatures are observed. Anomalies in co- 
inciding with specific heat anomalies are frequently ob- 
served in magnetic materials (see e.g. Refs. I50ll5lll52l ) . 
For an accurate phase-diagram determination to be de- 
scribed in the next sections, criteria for determining tran- 
sition temperatures are needed. In Fig. [3j we present 
temperature-dependent specific heat (a) between 110 and 
135 K, compared to ^^P- (b) and ^ for both current 
directions (c). There is a good agreement between the 
specific heat and (b), even though, strictly, the 

theories in Refs. [HI and are only good for paramag- 
netic to antiferromagnetic transitions. Clear anomalies, 
with (besides the fact that they point downward) similar 
shapes as the peak in specific heat, are also observed in 
^ (c) for both current-flow directions. 

In Fig. [31 we also present temperature-dependent spe- 
cific heat (d) between 195 and 230 K, compared to 
for both ac-field directions (e) and ^ for both current 
directions (f). Aside from the sharp downturn above ~ 
224 K, (see Fig. HJd), which we ascribed to an impurity 
of CeFe2, there is a good agreement between the specific 
heat (d) and ^^P- for both H _L c-axis, and reason- 
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able agreement for H j j c-axis(e). Such an agreement 
is expected for simple antiferromagnets as well as for 
less-simple antiferromagnets, see Refs.[53|and[54], respec- 
tively. As expected from the absence of a clear anomaly 
in the individual resistivity curves for both j / / c-axis 
and j _L c-axis, see Fig. [2fc, a clear anomaly is not ob- 
served in ^ (f) for either current-flow direction. 

Ideally, phase transitions are determined by specific 
heat. These measurements, however, are time-consuming 
and may be hindered by torque, and, as will be described 
below, by sample fragility. Therefore, for the determi- 
nation of phase transition fields and temperatures, we 
made use of magnetization and electrical resistivity data. 
Above, we found that Ce2Fei7 anomalies in the specific 
j p are also found in susceptibility via Fisher's re- 

2 dty g^ and in the resistivity via Below, we 
will make use of a modified version of Fisher's relation, 
d'dr^ ' an< ^ ver ify some of the results with from re- 
sistivity measurements, as well as with transition fields 
obtained from anomalies in obtained from magneti- 
zation isotherms^. 



V. PURE CE 2 FE 17 - RESULTS FOR H _L C-AXIS 



heat C r 
latior 



The plate-like samples, though, with their surface nor- 
mal parallel to the c-axis, were easy to (re)align, within 2- 
3°, and easier to clamp than the blocky samples. Because 
of their relative ease-of-use, we chose the plate-like sam- 
ples for the detailed magnetization measurements, and 
developed a protocol, see below, that helps circumvent 
their fragility. For determining hysteresis behavior, we 
used blocky samples. 




A. Magnetization 



2 3 4 5 
^ n H (T) 



We did not observe any noticeable differences in mag- 
netization for fields applied in different directions within 
the plane _!_ c-axis, therefore, they will be presented as 
measured with H _L c-axis. Before we discuss these ex- 
periments, we address the fragility of the compound. It 
has been reported^ that at ~ 5 K the crystallographic 
unit cell volume is larger in 5 T than in zero field, by 
~ 0.7 %. At 1.8 K, a hysteretic metamagnetic transi- 
tion, see below, takes place, presumably from a low-field, 
small-volume phase to a high-field larger-volume phase. 
In a magnetization experiment, at 1.8 K, on a plate-like 
sample, the sample shattered while the magnetic field (H 
_L c-axis) decreased through this transition. 

Although the data from the two types of crystal over- 
lap, we found plate-like samples to behave differently 
from blocky samples. On blocky samples, we have been 
able to measure full-hysteresis magnetization isotherms 
at temperatures between 1.8 K and 300 K. Visual inspec- 
tion after the measurements revealed no obvious damage. 
We do not know why plate-like samples appear more 
prone to shattering than the blocky samples. We may 
note, however, that the plate-like samples have a much 
larger surface-to- volume ratio than the blocky samples, 
and the relative change of the a and c lattice parameters 
reportedly^ is about equal. These effects can result in a 
relatively much larger relative change in surface area for 
the plate-like samples than for the blocky samples, when 
passing through the metamagnetic transition. Moreover, 
blocky samples may have more inclusions, which may 
prevent shattering. 



FIG. 4: Magnetization at various temperatures for Ce2Fei7 
with the magnetic field noH _L c-axis, measured with increas- 
ing and decreasing fields with steps of 0.1 T. For clarity, the 
curves have been offset with respect to one another. The black 
arrows indicate maxima in for increasing fields. 



Fig. [4] shows magnetization at various temperatures, 
measured on a blocky sample with both increasing and 
decreasing field strengths up to 7 T, applied _L the c- 
axis, with a step of 0.1 T between measurements. With 
increasing fields, the magnetization at 1.8 K is small and 
increases only slightly with increasing fields up to ~ 1 T. 
Then, around ~ 1 T, the magnetization jumps to sat- 
uration and is nearly field independent, reaching ~ 27 
/ie/f-u. at 7 T. This jump in magnetization shows a 
hysteresis of about 0.2 T. The magnetization at 50 and 
100 K also shows hysteretic jumps to saturation, at ~1.5 
and ~2.1 T, respectively. The width of the hysteresis 
becomes smaller with increasing temperature. Besides 
the hysteretic metamagnetic transition in high fields, the 
curves measured at 50 and 100 K both also show a step- 
like feature near 1 T with no detectable hysteresis, in- 
dicative of another metamagnetic transition. 

At temperatures above ~ 125 K, the lower zero-field 
transition temperature described in the previous section, 
magnetization isotherms are qualitatively different. At 
140 K, besides a step in magnetization that occurs here 
at ~ 1.2 T, the magnetization follows an s-shape increase 
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with increasing fields with a bending point near 2.6 T, 
up to ^3.3 T, where the magnetization starts to satu- 
rate. Hysteresis is not observed. At 200 K, still in the 
magnetically ordered phase, a weak feature is observed 
near ~ 0.9 T, and a change in slope may be discerned 
near ~ 2.7 T, indicating that the magnetization starts to 
saturate above this field strength. 
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FIG. 5: Magnetization as a function of (decreasing) temper- 
ature for fields (fioH _L c-axis) between 0.5 and 4.5 T. The 
arrows indicate maxima in 
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To avoid shattering of the plate-like samples, and still 
measure magnetization at any field/temperature, we de- 
veloped specific field and temperature protocols strate- 
gies. As indicated above, the sample that shattered, 
did so when it reentered the low-temperature low-field 
phase through a hysteretic phase transition, such as 
those shown in Fig. [4] Therefore, our measurements 
were performed so as to avoid such a transition. For 
measurements of temperature-dependent magnetization, 
with temperatures reaching below ~ 100 K, we used the 
following protocol. In fields above ~ 1 T, measurements 
were performed upon cooling from temperatures higher 
than the ordering temperature (~ 210 K). At low tem- 
peratures, the field was increased to above ~4T, which, 
as we will see below, is higher than any observed phase- 
transition field. After that, the sample was warmed to 
above 210 K, where the field was removed. Since the 
lowest field at which a hysteretic phase transition takes 
place is ~ IT, measurement in applied fields lower than 
this could be performed upon cooling as well as upon 
heating. 

Fig. [5] shows the field evolution of temperature de- 
pendent magnetization, measured on a plate-like sam- 
ple, at temperatures between 240 K and 80 K, using 
this protocol. In a field of 4.25 T, the magnetization 



increases monotonically with decreasing temperatures, 
without any clear anomalies. In lower fields, multiple 
anomalies occur, associated with magnetic phase tran- 
sitions. Based upon the results of the previous section, 
where it was found that the peaks in the zero-field 
appear very similar to the peaks in specific heat, we 
assume that peaks in 91 -q^ delineate magnetic phase 

boundaries. Anomalies determined from peaks in , 
are indicated by means of arrows in Fig. [5] Traces of 
are shown in Fig. [5b, as well as the positions of 
the peaks and their width, determined from half width 
at half maximum on the steepest flank. 

Detailed measurements of field-dependent magnetiza- 
tion at fixed temperatures were performed using a similar 
protocol as above. For each temperature, 5 K, and 10 K 
up to 220 K with 10 K steps, the sample was cooled in 
zero field, and the field was increased in steps of 0.025 T 
up to the highest measuring field. This field was chosen 
such that the magnetization appeared clearly saturated. 
After that, the field was increased up to above 4.5 T, 
without changing the temperature, and then the tem- 
perature was increased to above 250 K, at which temper- 
ature the field was removed. The results of these mea- 
surements are shown in FigfBJa, as a contour plot, that 
displays the magnetization as a function of applied field 
and temperature. 

As an ofte n 51 i 55 i 56 applied criterion for transition 
fields, according to mean-field theory (See e.g. [55II56T ). 
peaks in the differential susceptibility determined 



from magnetization isotherms, can be taken. Fig. 

3M 
dH 



shows several traces of The peaks, and their width, 



determined from the half width at half maximum on the 
steepest flank, are also indicated in Fig. [Sh- The differ- 
ences between the peak positions indicated here and in 
Fig. [4] can be ascribed to differences in demagnetizing 
fields, due to the different shapes of the samples. 

Displayed together on top of the magnetization 
isotherms, in Fig. [SJd, the positions and the widths of 
the peaks obtained from both 1^ and from , ap- 

pear to overlap. We therefore conclude that they indicate 
magnetic phase boundary lines, and thus outline a mag- 
netic phase diagram. The phase boundary line that is 
associated with clear hysteresis is drawn in grey. For the 
other phase boundary lines, hysteresis was not observed. 
The labels will be discussed below. 

To determine details of the apparent crossing of phase 
boundary lines at about 120 K and in about 1 T, we 
performed a more detailed experiment. Temperature- 
dependent magnetization was measured upon cooling at 
1 K steps between 140 K and 100 K. Fields ranged be- 
tween 1.5 T and 0.925 T at 0.025 T steps. The results of 
the measurements are displayed in Fig. GJl, as a contour 
plot, that displays the magnetization as a function of ap- 
plied field and temperature. In Fig. [TJa, we display the 
derivative 8 ^ T ^ as a contour plot. The peak positions 
and their widths, obtained from this experiment are in- 
cluded in Fig. [7^., as well as peak positions and widths 
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FIG. 7: Detail of H-T phase diagram for H _L c-axis. a Con- 
tour plot of magnetization M(H,T) with peaks determined 
from 8 ^ T ^ (horizontal bars) and peaks determined from 4^ 



(vertical bars) plotted on top. b Contour plot of 
termined from a. 



B. Resistivity 
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FIG. 6: Magnetic phase diagram for for Ce2Fei7 obtained 
with the applied magnetic field H _L c-axis. a Traces of 
obtained from field-dependent magnetization. Peaks and 
width of peaks are indicated, b Contour plot of magnetiza- 
tion M(H, T) with peaks from a and c plotted on top. Labels 
are explained in the text. Dark grey indicates where no data 
was taken, c Traces of 8 ^q"^ T ' > obtained from temperature- 
dependent magnetization. Peaks and width of peaks are in- 
dicated. 



obtained from the magnetization isotherms in Fig. [6}d. 
From this it is clear that, within experimental resolution, 
there is a crossing of phase boundary lines at ~ 122 K in 
- 1.2 T. 



We determined temperature-dependent resistivity for 
the current j flowing / / c-axis as well as for j _L c-axis 
in various applied fields H _L c-axis up to 5.5 T. Af- 
ter the zero-field measurements, we performed measure- 
ments in applied fields. With the fragility of the material 
in mind, we measured, starting at 300 K, with both de- 
creasing and increasing temperatures. In this way, we 
verified that the samples remained intact, which was ob- 
vious if, at a given temperature, the resistance measured 
upon heating was found to be the same as the resistance 
measured upon cooling. The measurements with j / / c- 
axis were performed at gradually lower fields, starting at 
5.5 T. Down to 2.25 T, there was no significant change 
in the resistance determined with increasing or with de- 
creasing temperatures. At 2 T, however, upon cooling 
below ~ 86 K, the sample resistance became unstable: 
it started to fluctuate, and the resistance upon warming 
was higher than upon cooling. Note that at this tem- 
perature and field the phase diagram, Fig. O shows a 
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FIG. 8: a Temperature-dependent normalized resistivity FIG. 9: a Traces of for j // c-axis, plot on top of magnetic 
measured with j // c-axis in different H _L c-axis up to 5.5 T. phase diagram for H _L c-axis. b as in a, but for j _L c-axis. 
b as in a, but for j J_ c-axis. 



hysteretic transition. Although visual inspection did not 
reveal any obvious damage, the sample was no longer us- 
able. For the magnetization measurements, we were able 
to prevent the samples from shattering by crossing the 
hysteretic phase boundary only with increasing field or 
decreasing temperature. Apparently, for the resistivity 
measurements, such a hysteretic phase transition has to 
be completely avoided. For this reason the lowest tem- 
perature for the measurement in 1.5 T with j _L c-axis 
was 80 K. Besides these, for the other applied- field val- 
ues, measurements were performed every 1 K between 
300 K and 2 K. 

The results are shown in Fig. [5^ for j / / c-axis, and 
in Fig. [SJa for j J_ c-axis. For completeness, both Fig. [5^ 
and b include the zero-field results from FigO Figs. [S^ 
and b shows traces of ^ plot over the phase diagram 
from Fig. [SI Also included in Fig.[S]are bars that indicate 
peak positions and widths in ^ . Note that the peaks in 
^ coincide well with the previously determined phase 
diagram. 

The anomalous increase in resistivity with decreasing 
temperature, effects we ascribed above to the opening 
of an anisotropic superzone gap, observed in zero field 



for both current directions, occurs also in applied fields 
for both current directions. We observed it for fields of 
2 T (j j! c-axis), 1.5 T, 0.5 T (j 1 c-axis). Note that in 
zero field and in 0.5 T, the crossed boundary lies between 
phases marked as AF1 and AF2, in Figs.[S]and[71 whereas 
for 1.5 T and 2 T, the crossed boundary lies between the 
phases marked as SF and SF2. 

Generally, for both employed current directions, the 
high-field resistivity at low temperatures is low, whereas 
for low fields it is high. The residual resistance ratio 
in 5.5 T, about 14 and 17, for j // c-axis and j _L 
c-axis, respectively, is much higher than in zero field. 
Since, as already indicated above, in Sec. IIV1 the resis- 
tivity for j _L c-axis is hardly affected by superzone-gap 
effects, whereas its low-temperature resistivity is much 
lower in 5.5 T than in zero-field, the superzone gap does 
not determine the large magnetoresistance, as was previ- 
ously noted^. Assuming the resistivity is determined as 
p = po + piatticc + Pmag, and that the residual resistiv- 
ity, po, and the lattice contribution piattice are the same 
for the low-temperature AF2 or P phase, we consider 
the large magnetoresistance mainly caused by a strong 
change in magnetic scattering /? m ag- 

For j _L c-axis the anomalies associated with phase 
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transitions are weaker than for j / / c-axis. The ordering 
transition between the paramagnetic phase P and AF1 
(for 0.5 T) or between P and SF3 (for 1.5 and 3 T) is 
not observed as an anomaly for j _L c-axis, nor is the 
transition between SF3 and SF1 (for 1.5 T). For j 1 c- 
axis, the transition between SF3 and P produces a clear 
anomaly. For j j j c-axis, weak anomalies mark the phase 
boundaries between P and SF3 (at 2 T, 2.5 T and 3 T). 
For the purpose of this paper, we do not consider the 
'break' in resistivity around 175 K, for j / / c-axis in 
5.5 T, resulting in a very broad 'peak' in due to a 
phase transition. 

Although the presented dataset is limited, we have 
been able to measure resistivity in every (H, T) phase 
in the phase diagram of Fig. [6l Our results indicate 
that the previously observed large magnetoresistance~ 
is mainly due to a large (magnetic) scattering of conduc- 
tion electrons in the low-temperature AF1 phase, which 
is strongly reduced in the field-induced P state. Here, 
the superzone gap is associated with the phase-boundary 
line that separates AF1 and SF1 (above 124 K), from 
AF2 and SF2. 



VI. PURE CE 2 FEi7- MAGNETIZATION 
RESULTS FOR H // C-AXIS 
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FIG. 10: Magnetization at various temperatures for Ce2Fei7 
with the magnetic field /j,qH // c-axis, measured with increas- 
ing and decreasing fields with steps of 0.1 T. For clarity, the 
curves have been offset with respect to one another. The black 
arrows indicate maxima in for increasing fields. 



In this section we discuss magnetization for the field 
applied // c-axis, resulting in an H-T magnetic phase di- 
agram for this applied-field direction. As in Sec. IV Al we 
used blocky samples for determining hysteresis behavior, 
and we used plate-like samples for detailed magnetization 
measurements. 

Fig. [TU] shows magnetization at various temperatures 
up to 145 K, measured on a blocky sample with both 
increasing and decreasing field strengths up to 5.5 T, ap- 
plied // c-axis, with a step of 0.1 T between measure- 
ments. With increasing fields, the magnetization at 2 K 
is small and increases slightly with increasing fields up to 
~ 3.2 T. Then, around ~ 3.2 T, a metamagnetic tran- 
sition takes place, followed by a second metamagnetic 
transition at ~ 3.9 T . Both these transitions show hys- 
teresis. At 5.5 T the magnetization is not saturated, and 
it reaches a value of ~ 26 /ie/f-u. lower than the magne- 
tization observed at 1.5 T for H _L c-axis, in Fig.0J The 
lower-field metamagnetic transition occurs at tempera- 
tures below ~ 125 K, at which a zero-field transition oc- 
curs. The (increasing) field, at which this metamagnetic 
transition occurs, decreases slightly up to 80 K, and at 
the same time the hysteresis gradually becomes smaller. 
The field at which the higher-field transition occurs, in- 
creases with increasing temperature, at least up to 80 K. 
The hysteresis associated with this transition is observed 
at 20 K, but not at 40 K, nor at higher temperatures. At 
145 K, the magnetization increases smoothly, and with- 
out anomalies, until it approaches saturation near 5.5 T. 

Fig. ITT1 shows the field evolution of temperature depen- 
dent magnetization, measured on a plate-like sample, at 
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FIG. 11: Magnetization as a function of (decreasing) temper- 
ature for fields (/lo-ff // c-axis) between 0.5 and 7 T. The 
arrows indicate maxima in 
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temperatures between 240 K and 60 K. A similar strat- 
egy as described above in Sec. IV Al was used. We mea- 
sured upon cooling, and at 60 K, we increased the field 
up to 7 T, and heated the sample to ~ 240 K, before we 
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removed the field. In a field of 7 T, the magnetization in- 
creases uniformly with decreasing temperatures, without 
clear anomalies. In lower fields multiple anomalies oc- 
cur, associated with magnetic phase transitions. Based 
upon the results of the previous sections, we assume that 
peaks in a ^ T - > in applied fields denote magnetic phase 

boundaries. Anomalies determined from peaks in a ^ T ' 1 , 
are indicated by means of arrows in Fig. 1111 Traces of 
a ^g I j^ > are shown in Fig. [T2"b . as well as the positions of 
the peaks and their width, determined from half width 
at half maximum on the steepest flank. 

Detailed measurements of field-dependent magnetiza- 
tion at fixed temperatures were performed using a similar 
protocol as above. For each temperature, 5 K, and 10 K 
up to 220 K with a 20 K step, the sample was cooled in 
zero field, and the field was increased in steps of 0.05 T 
up to 7 T. Then the temperature was increased to above 
250 K, at which temperature the field was removed. The 
results of these measurements are shown in Fig |12b . as a 
contour plot, that displays the magnetization as a func- 
tion of applied field and temperature. 

The anomalies observed in these magnetization 
isotherms produced peaks in the derivatives Fig.[T2"b 
shows several traces of The peaks, and their width, 
determined from the half width at half maximum on the 
steepest flank, are also indicated in Fig. [T2k The differ- 
ences between the peak positions indicated here and in 
Fig. [10] can be ascribed to differences in demagnetizing 
fields, due to the different shapes of the samples. 

As for Fig. [6j in Fig. [12b the positions and the widths 
of the peaks obtained from both and from a ^ T - 1 , 
overlap. We therefore conclude that they indicate mag- 
netic phase boundary lines, and thus outline a magnetic 
phase diagram. Phase boundary lines associated with 
hysteresis are indicated with gray lines. For the other 
phase boundary lines, hysteresis was not observed. The 
labels will be discussed below. 

Fig. [T3] shows details of the apparent splitting of 
phase boundary lines at about 100 K and in about 
3.2 T. Temperature-dependent magnetization was mea- 
sured upon cooling at 1 K steps between 140 K and 60 K. 
Fields ranged between 3.5 T and 3.0 T at irregular inter- 
vals. The results of the measurements are displayed in 
Fig. 113b. as a contour plot, that displays the magnetiza- 
tion as a function of applied field and temperature. The 
contour plot shows field steps of 0.025 T, values for mag- 
netization in fields between the measurement fields were 
obtained by interpolation. In Fig. 113b . we display traces 
of the derivative ■ The peak positions and their 

widths, obtained from this experiment are included in 
Fig.[T3"b. Fig.[T3"b also contains these peak positions and 
their widths, as well as widths obtained from the mag- 
netization isotherms in Fig. 112b . From this it is clear, 
that within experimental resolution, there is a splitting 
of phase boundary lines at ~ 100 K in ~ 3.15 T. 
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FIG. 12: Magnetic phase diagram for for Ce2Fei7 obtained 
with the applied magnetic field H // c-axis. a Traces of 
obtained from field-dependent magnetization. Peaks and 
width of peaks are indicated, b Contour plot of magnetization 
M (H, T) with peaks from a and c plotted on top. Labels 
are explained in the text, c Traces of 9( -gy T - > obtained from 
temperature-dependent magnetization. Peaks and width of 
peaks are indicated. 



VII. PURE CE2FE17- ANGLE-DEPENDENT 
MAGNETIZATION 



Magnetization for fields applied at different angles 9 
with the c-axis, was measured on a plate-like sample 
(mass < 0.5 mg), mounted on a rotator. We measured 
over a range of 130° with 5° steps, overshooting both H 
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FIG. 13: Detail of H-T phase diagram for for H / / c-axis. a 
Traces of a< -Qy T - ) obtained from temperature-dependent mag- 
netization. Peaks and width of peaks are indicated as horizon- 
tal bars, b Contour plot of magnetization M (H, T), obtained 
via interpolation, with peaks determined from 9 (hori- 
zontal bars) and peaks determined from (vertical bars) 
plotted on top. 



FIG. 14: Evolution of magnetization isotherms, at 100 K (a, 
and at 150 K (b), measured at angles with the c-axis, with a 
15° step, from / / c-axis, to _L c-axis. For clarity the curves 
are offset with respect to one another. The arrows indicate 
extrema in ^M. 



II c-axis and H _L c-axis. As above, to preserve the sam- 
ple, we used a protocol. In zero field, and at ~ 240 K, 9 
was set. Then the sample was cooled to the measurement 
temperature, at which the magnetization was measured 
with increasing fields, every 0.05 T up to 7 T. After that, 
the sample was heated up, in 7 T, to ~ 240 K, at which 
temperature the field was removed, and a new 9 was set. 

The evolution of the magnetization, as 9 is varied, for 
both T = 100 K and T = 150 K, is repectively shown in 
Fig. 114b and in Fig. 114b . Arrows indicate extrema found 
in Qfjj-, note that some extrema are more clearly visible 
than others (cf. Figs.fflandUD]). For H // c-axis, at 150K, 
th magnetization increases smoothly, up to about 5.5 T, 
where the magnetization starts to saturate. No s-shapes 
marking phase transitions were observed for H / / c-axis. 
As 9 increases, near 1-2 T an s-shape appears, that grad- 
ually develops into the step-like transition observed for 
H _L c-axis. The transition that marks the field where 
the magnetization starts to saturate, gradually decreases 



with increasing 9. Furthermore, with 9 increasing above 
~ 30 °, an s-shape appears at fields slightly lower than 
the saturation field. 

Also at 100 K, with increasing 9, the field at which 
the magnetization starts to saturate decreases gradu- 
ally. However, here it develops from a smooth, slightly 
s-shaped, non-hysteretic, transition into a very sharp and 
hysteretic transition, see Figs. [10] and [4] respectively. 
With increasing 8, the sharp (and hysteretic) transition 
observed for H // c-axis near 3 T shifts to a slightly 
lower field, and for 8 > 45°, appears to merge with the 
jump-like transition to saturation observed for H _L c- 
axis. Furthermore, the step-like transition near 1.5 T, 
observed for H _L c-axis develops from an s-shape that 
first appears for 8 > 30°. 

Contour maps of M(8,H) measured at 100 (a) and 
150 K (b), which is below and above the 125 K zero-field 
transition, respectively, are displayed in Fig. 1151 Vertical 
bars denote peaks in I ^C a nd their widths obtained as 
above in Sees. IV Al and IVIl Note that, especially for the 
less-visible transitions, the peak widths may be slightly 
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FIG. 15: Magnetic phase diagrams determined at 100 K and 
150 K for fields applied at various angles 6 with the c-axis. a 
Contour plot of magnetization M(H,8) at 100 K with peaks 
determined from (vertical bars) plotted on top. Labels 
are explained in the text, b Contour plot of magnetization 
M (H,9) at 150 K with peaks determined from ^1 (vertical 
bars) plotted on top. Labels are explained in the text. 




underestimated (cf. Figs. [6l and [T2|) . For Fig. [T5h. mea- 
sured at 100 K, for H _L c-axis, transitions were found 
near ~ 1 T and near ~ 2 T, just as in Fig. [51 At the same 
temperature, for H / / c-axis, transitions were found near 
- 3 T and near ~ 5.7 T, just as in Fig. [H Fig. HBJo, 
measured at 150 K, compares similarly to Figs.[S]and [T2] 
The labels are discussed in Sec. IIXI 



VIII. TA-DOPED CE2FE17- MAGNETIZATION 
RESULTS FOR H 1 C-AXIS 



FIG. 16: a Magnetization at various temperatures for Ta- 
doped Ce2Fei7 with the magnetic field hqH _L c-axis, mea- 
sured with increasing fields with steps of 0.1 T. For clarity, 
the curves have been offset with respect to one another. The 
black arrows indicate maxima in . b Contour plot of mag- 



netization M (H, T) with peaks from ^£ and 



plotted 



8T 

on top. Labels are explained in the text, c Temperature- 
dependent magnetization in various fields up to 5 T. Arrows 
indicate phase-transition temperatures. 



In order to compare undoped Ce2Fei7 to Ta-doped 
Ce2Fei7, we used the Ta-doped Ce2Fei7 crystal with 
nominally 0.5 % Ta, that was also used for the EPMA ex- 
periment. We chose to measure magnetization with the 
field H _L c-axis, because, as will become clear below, 
this is the most salient applied-field direction. Fig. [T|5] 
shows a summary of these measurements. Fig. 116b shows 
temperature-dependent magnetization curves, measured 



in 0.05, 0.5, 1, 1.5, 2, 3, 4, and 5 T, at (decreasing) tem- 
peratures between 250 K and 5 K. In curves of 8 ^ T - > 
(not shown), extrema were found, which are indicated by 
arrows in Fig. 116b . In fields of 3 T and higher, no sharp 
anomalies were observed, similar to undoped Ce2Fei7 for 
H 1 c-axis in 4.25 T (see Fig.©. In the lower field of 2 T, 
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anomalies appear, marking phase transitions. In contrast 
to undoped Ce2Fei7, a large magnetization is found be- 
low ~ 75 K, in all applied fields of 0.05 T and higher, 
indicating that our sample of Ta-doped Ce2Fei7 is ferro- 
magnetic at low temperatures. The multiple anomalies 
in fields higher than 0.05 T, especially in 1 T, hint at a 
complex phase diagram. The low-field 0.05 T results are 
consistent with two ordering transitions, antiferromag- 
netic at Tn = 214(1), and ferromagnetic at To = 75(5) K. 
Results of (increasing) field-dependent magnetization at 
various temperatures between 10 K and 200 K are dis- 
played in Fig. [TOb . At 10 K and at 50 K the magnetiza- 
tion shows ferromagnetic behavior, and no metamagnetic 
transitions are observed. The saturation magnetization 
of ~ 31 /ie/f-u. at 10 K is in good agreement with results 
published for Ta-crucible grown Ce2Fei?22. At 100 K, 
150 K, and 200 K, metamagnetic transitions were ob- 
served as maxima in at positions which are indicated 
by arrows. At 100 K, only one anomaly was observed, 
whereas at 150 and 200 K two anomalies were observed, 
indicating the existence of a metamagnetic phase at these 
temperatures, and consistent with the multiple anomalies 
found in temperature-dependent magnetization. Fig. 116b 
shows a contour map of M(H,T), and phase boundaries 
obtained from (vertical bars) and ^g^P (horizontal 
bars). The labels are analogous to those given for pure 
Ce2Fei7, but here is only one zero-field AF phase. 



IX. DISCUSSION 

The detailed study of the magnetic phase diagrams of 
undoped and of Ta-doped Ce2Fei7 provides insight into 
the competition between ferro- and antiferromagnetic in- 
teractions in these compounds. In a magnetic phase dia- 
gram, one can find a phase-boundary line starting, in zero 
field, at the Neel temperature. It separates the paramag- 
netic phase from the magnetically ordered phase. Gener- 
ally, its critical magnetic field strength increases with de- 
creasing temperature 5 - 5 -. However, for Ta-doped Ce2Fei7, 
the ground state is not antiferromagnetic but ferromag- 
netic. In such a case, the antiferromagnetic transition 
fields may form a 'dome' in H-T, as predicted in Ref . [57l. 
Such a dome is indeed observed in the H-T phase dia- 
gram of Ta-doped Ce2Fei7, which in zero fields orders AF 
at the Neel temperature Tn ~ 214 K, and ferromagnetic 
(F) at Tb~ 75 K. Furthermore, besides the AF dome, 
there is an additional dome, enclosing a metamagnetic 
spin-flop-like phase, SF. This dome has a maximum field 
of - 2 T near 170 K. 

Undoped Ce 2 Fei 7 also order AF, at T N ~ 208 K. 
Above T t ~ 124 K, also the H-T phase diagram for H 
_L c-axis Fig O shows an AF dome. Besides this, at 
there are two more domes, enclosing metamagnetic spin- 
flop-like phases, SF1 and SF3. The dome enveloping 
SF3 has a maximum field of ~ 4 T at ~ 170 K. How- 
ever, the boundary of this high-temperature SF3 dome 
does not go down to zero field at some lower tempera- 



ture, instead it converges with another, low-temperature, 
dome, which is connected to T t at ~ 124 K, which is, as 
reported elsewhere in this volume^, associated with a 
crystallographic phase transition. This low-temperature 
dome also encloses a metamagnetic spin-flop-like phase, 
labeled SF2. The result is that undoped Ce2Fei7 remains 
antiferromagnetic below T t (AF2). 

As demonstrated in Fig. there is a crossing of phase 
boundary lines at ^122 K in ^1.2 T. According to Lan- 
dau^, for continuous phase transitions such crossings are 
only possible if the crossing phase boundary lines rep- 
resent different order parameters. Although from mean- 
field descriptions it may be concluded that spin-flop tran- 
sitions are first order—, we observed no hysteresis in the 
measurements of magnetization around the phase bound- 
ary line separating the low-field phases AF1 and AF2 
from the spin-flop-like phases SF1 and SF2, respectively. 
The phase boundary that starts at T t in zero fields is 
characterized by a superzone-gap-like feature in electri- 
cal resistivity. The superzone-gap-like behavior was ob- 
served for the phase transition from AF1 to AF2, and 
between SF1 and SF2. We therefore assume that the 
field-dependent superzone-gap-like behavior is connected 
to the observed zero-field crystallographic phase transi- 
tion^. 

The low-temperature dome that starts at T t in zero 
field ends in about 1 T at K. The boundaries of the 
low-temperature dome indicated by a grey line in Fig. [5] 
are due to hysteretic phase transitions. The boundaries 
of the high-temperature dome do not show hysteresis. 
The high-temperature dome is associated with an anti- 
ferromagnetic phase transition, and the low-temperature 
dome with a crystallographic phase transition. Accord- 
ing to Landau 5 -^, a convergence of two phase-boundary 
lines of continuous phase transitions may produce a 
line of first-order phase transitions, under certain con- 
ditions 59 . In this picture, the grey line of hysteretic first- 
order phase transitions constitutes a combined crystallo- 
graphic and antiferromagnetic phase transition. 

Now we turn to the H- T phase diagram H / / c-axis 
for undoped Ce2Fei7, Fig.[T2J At the Neel temperature, 
a phase boundary line between AF1 and P starts. It 
is connected to a phase boundary line enveloping SF4, 
and a phase boundary line between SF4 and AF1. The 
boundary between SF4 and P becomes hysteretic below 
~ 20 K. If the s-shape magnetization just below 140 K 
is due to a continuous phase transition, then this phase 
boundary line may have a tricritica l 55 ! 56 point near 20 K, 
see Fig [TO] 

Connected to T t is a phase-boundary line, between 
AF1 and AF2, that is associated with a crystallographic 
phase transition^, this line is connected to the hysteretic 
phase-boundary line between AF2 and SF4, and to the 
phase-boundary line between AF1 and SF4. Obviously, 
the hysteretic phase boundary line between AF2 and SF4 
represents first-order phase transitions. Then, follow- 
ing thermodynamic argument s 58 ! 59 , the phase-boundary 
line between AF1 and AF2, and the phase-boundary line 
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between AF1 and SF4 may represent continuous phase 
transitions at this point. 

Comparing the magnetization behavior near the T t - 
line in fields below ~2.5 T, it appears that the magneti- 
zation at lower temperatures is lower than at higher tem- 
peratures. Thus the crystallographic phase transition 32 
affects the magnetic properties for this applied field direc- 
tion. The situation is different for higher fields, because 
at a given field above ~3.5 T, the magnetization in the 
low-temperature phase SF4 is higher than in the neigh- 
boring phase AF1. This suggests that the phase bound- 
ary line between AF1 and SF4 is associated with the 
occurrence of a low-temperature field-stabilized phase. 

At 150 K, well above Tt, the antiferromagnetic phase 
transition is observed for both H // c-axis and for H _L 
c-axis near 3.5 T and 6 T respectively, whereas spin- 
flop-like phase-transitions are only observed for H // 
c-axis. The angle-dependent magnetization experiment, 
Fig. 115b shows how the antiferromagnetic transition and 
the transitions between AF1 and SF1, and between SF1 
and SF3, respectively, evolve, as the angle between the 
applied field and the c-axis is changed. In simple cases^, 
spin-flop transitions are observed for weakly anisotropic 
antiferromagnets, when the field is applied parallel to 
the ordered moment direction. In this picture, the angle- 
dependent results measured at 150 K indicate that the 
moments in Ce2Fei7 order in the plane perpendicular 
to the c-axis. Moreover, the angle-dependent behavior 
of the antiferromagnetic phase-transition field, which is 
higher for H // c-axis than for H _L c-axis, is also indica- 
tive of a magnetic anisotropy in favor of the plane _L c- 
axis. Also at 100 K, well below T t , the angle-dependence 
of the antiferromagnetic phase-transition field, higher for 
H 1 1 c-axis than for H _L c-axis, indicates a magnetic 
anisotropy in favor of the plane _L c-axis. The pres- 
ence of the non-hysteretic spin-flop-like transition be- 
tween AF2 and SF2, for H J_ c-axis, and its absence 
for H If c-axis, is also an indication that at 100 K 
the ordered magnetic moments lie in the plane _L c- 
axis. The angle-dependent magnetization furthermore 
indicates that SF4, observed for H // c-axis, is not con- 
nected to any phase that can be observed for H _L c-axis. 
Therefore the angle-dependent results are consistent with 
the interpretations that the hysteretic AF2-SF4 phase- 
boundary line includes the crystallographic phase tran- 
sition, and that the phase-boundary line enveloping SF4 
is related to the magnetic anisotropy, as also observed at 
150 K. 

We think the H-T phase diagrams of Ta-doped 
Ce2Fei7 and of undoped Ce2Fei7 are manifestations 
of competing antiferromagnetic and ferromagnetic ex- 
change interactions. The occurrence of magnetic phase 
transitions in itinerant-electron systems with one type 
of ferromagnetism and one type of antiferromagnetism 
in competition, has been treated theoretically by Moriya 
and Usami^i. They constructed possible phase diagrams, 
one of which is very similar to the Ta-doped Ce2Fei7 for 
H _L c-axis. This predicted phase diagram has a low- 



temperature ferromagnetic state, an antiferromagnetic 
phase at intermediate temperatures, and a paramagnetic 
state at high temperatures. In non-zero applied fields, 
the antiferromagnetic phase is enveloped by a dome. 

Although similar, the phase diagram for Ta-doped 
Ce2Fei7, Fig.Ql)]is more complex than the one described 
by Moriya and Usami, because here the AF dome shows 
another dome on top of it, indicated as SF. In the work 
of Moriya and Usami, calculations are made for two 
competing and interacting magnetic modes, a uniform 
magnetization Mq, and a staggered magnetization Mq, 
for the competing ferromagnetic state with wavevector 
and antiferromagnetic state with wave vector Q, re- 
spectively. The free energy was expanded up to the 
fourth order in these magnetic modes, without consid- 
ering magnetic anisotropy. It therefore seems that the 
model can be adjusted by including higher-order terms 
in magnetization, by including magnetic anisotropy or 
by introducing more than two competing types of ex- 
change interactions. Although the, isotropic, Moriya- 
Usami model has not been extended with an in-plane 
easy-plane anisotropy, no in-plane anisotropy has been 
observed for Ce2Fei7, and easy-plane systems, with mo- 
ments confined to the easy plane are expected to behave 
similarly to isotropic systems. Higher order terms in 
magnetization are not expected to qualitatively change 
the possible phase diagrams^. We think that an exten- 
sion of the theory by possibly including more than two 
competing and interacting magnetic modes may result in 
a phase diagram similar to the one we observed for Ta- 
doped Ce2Fei7. This seems plausible, in view of the com- 
plexity of the crystal structure, with 5 different crystallo- 
graphic sites. Moreover, electronic-structure calculations 
on related Y2Fei7 produced different exchange interac- 
tion parameters of both antiferromagnetic and ferromag- 
netic type, between these siteai^. For Al-doped Fe3Ga4, 
Duijn et ali^ described a similar H-T phase diagram in a 
similar manner. In such a picture, the phase SF is canted, 
with a non-zero magnetization, and, moreover, is char- 
acterized by a, possibly additional, wave vector different 
from the zero-field wave vector. This could possibly be 
verified by a neutron-diffraction experiment. 

In a similar way, the H-T phase diagram for undoped 
Ce2Fei7 can be evaluated. Here, at temperatures above 
Tt, there exists, besides the zero- field phase AF1, two ad- 
ditional ordered phases, SF1 and SF3. If the introduction 
of an additional competing magnetization for Ta-doped 
Ce2Fei7 holds, both SF1 and SF3 are canted, and even 
more modes Mq are relevant here. Within the frame- 
work of the Moriya-Usami theory, a small amount of (Ta) 
dopant influences the competition of Mq and relevant 
Mqs and the number of relevant Mq in the intermediate- 
temperature ordered state. As indicated in the introduc- 
tion, the electronic density of states is probably strongly 
peaked near the Fermi-surface, and therefore the elec- 
tronic and magnetic properties are likely strongly depen- 
dent on volume. Different intermediate-temperature be- 
havior of Ta-doped Ce2Fei7 and undoped Ce2Fei7 may 
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well be mostly related to differences in volume. 

X. CONCLUSIONS 

The magnetism of Ce2Fei7, earlier thought to be ab- 
normal, may be explained as being due to competing 
fcrro- and antiferromagnetism, and is well described by 
a Moriya-Usami-type phenomenological theory, which in 
its unmodified form yields magnetic phase diagrams for 
systems with a competition between one antiferromag- 
netic and one ferromagnetic mode. Magnetic phase di- 
agrams obtained for undoped Ce2Fei7 indicate that the 
magnetic order is determined by the competition between 
multiple antiferromagnetic modes and a ferromagnetic 
mode. This competition and the resulting magnetic or- 
der are sensitive to minor changes in the crystallographic 
unit cell. A minor amount of Ta preferentially substi- 
tutes for Fe on one of the 4 Fe crystallographic sites and 
enlarges the unit-cell volume. This Ta doping influences 
both the competition and the number of competing mag- 
netic modes, and results in a ferromagnetic ground state. 

Within the ordered state of undoped Ce2Fei7, as de- 
tailed elsewhere in this volume^, a doubling of the unit 
cell in the c-direction takes place, which is due to a minor 
atomic displacement. Also this crystallographic mod- 
ification has a distinct effect on both the competition 



and the number of competing magnetic modes. Here, 
the resulting ground state is antiferromagnetic, and it 
may thus be said that the crystallographic phase transi- 
tion stabilizes the antiferromagnetism. The sensitivity of 
the magnetism to Ta doping and to a minor displacive 
change in crystal structure is likely related to a strong 
sensitivity of electronic properties to small changes in 
(local atomic) volume, which in turn is consistent with a 
strongly peaked density of states near the Fermi level, as 
reported for the related compound Y 2 Fei 7 . 

In conclusion, the present work indicates that the same 
mechanisms that determine the magnetic properties of 
R-2Fei7 compounds, the presence of both ferro- and anti- 
ferromagnetic exchange interactions, and lattice anoma- 
lies, play a crucial role in determining the magnetism in 
Ce 2 Fei 7 . 

Acknowledgments 

The authors are indebted to S. L. Bud'ko, R. W. Mc- 
Callum, K. W. Dennis, J. Frederick, S. Jia, M. Angst, 
R. J. McQueeney, and A. I. Goldman. Ames Laboratory 
is operated for the U.S. Department of Energy by Iowa 
State University under Contract No. W-7405-ENG-82. 
This work was supported by the Director for Energy Re- 
search, Office of Basic Energy Sciences. 



* Email: yjanssen@ameslab.gov 

1 K. H. J. Buschow, Rep. Prog. Phys. 54, 1123 (1991). 

2 K. H. J. Buschow and J. S. van Wieringen, Phys. Stat. Sol. 
p. 231 (1970). 

3 K. H. J. Buschow, Rep. Prog. Phys. 40, 1179 (1977). 

4 K. H. J. Buschow, J. Less-Common Met. 11, 204 (1966). 

5 D. Givord and R. Lemaire, IEEE Trans. Magn. MAG-10, 
109 (1974). 

6 Y. Janssen, H. Fujii, T. Ekino, K. Izawa, T. Suzuki, T. Fu- 
jita, and F. R. de Boer, Phys. Rev. B 56, 13716 (1997). 

7 D. P. Kozlenko, V. P. Glazkov, B. N. Savenko, and V. I. 
Voronin, Eur. Phys. J. B 41, 445 (2004). 

8 J. Kamarad, O. Prokhnenko, K. Prokes, and Z. Arnold, J. 
Phys.: Condens. Matter 17, S3069 (2005). 

9 O. Prokhnenko, J. Kamarad, K. Prokes, Z. Arnold, and 
A. V. Andreev, Phys. Rev. Lett. 94, 107201 (2005). 

10 R. Coehoorn, Phys. Rev. B 39, 13072 (1989). 

11 D. Givord and R. Lemaire, Proceedings of ICM 197,3 
(Nauka, Moscow, 1974), p. 492. 

12 R. F. Sabiryanov and S. S. Jaswal, Phys. Rev. B 57, 7767 
(1998). 

13 L. M. Sandratskii, Adv. Phys. 47, 91 (1998). 

14 K. Knopfle, L. M. Sandratskii, and J. Kiibler, Phys. Rev. 
B 62, 5564 (2000). 

15 M. Marsman and J. Hafner, Phys. Rev. B 66, 224409 
(2002). 

16 S. Shallcross, A. E. Kissavos, S. Sharma, and V. Meded, 
Phys. Rev. B 73, 104443 (2006). 

17 R. Plumier and M. Sougi, Proceedings of ICM 1973 
(Nauka, Moscow, 1974), p. 487. 



18 K. Koyama, T. Goto, H. Fujii, N. Takeshita, N. Mori, 
H. Fukuda, and Y. Janssen, J. Phys. Soc. Jpn. 67, 1879 
(1998). 

19 K. Koyama, H. Fujii, T. Goto, H. Fukuda, and Y. Janssen, 
Physica B 294, 168 (2001). 

20 F. Ishikawa, T. Goto, and H. Fujii, Physica B 329, 329 

(2003) . 

21 Y. Makihara, H. Fujii, T. Fujiwara, K. Watanabe, K. Taka- 
hashi, K. K. Koyama, and M. Motokawa, Physica B 329, 
329 (2003). 

22 Y. Makihara, Y. Uwatoko, H. Matsuoka, M. Kosaka, 
H. Fukuda, and H. Fujii, J. Magn. Magn. Mater. 272, 551 

(2004) . 

23 O. Prokhnenko, C. Ritter, Z. Arnold, O. Isnard, J. Ka- 
marad, A. Pirogov, A. Teplykh, and A. Kuchin, J. Appl. 
Phys. 92, 385 (2002). 

24 P. A. Teplykh, A. Pirogov, A. Kuchin, and A. Teplykh, 
Phys. B 350, e99 (2004). 

25 D. P. Middleton, S. R. Mishra, G. J. Long, O. A. Pringle, 
Z. Hu, W. B. Yelon, F. Grandjean, and K. H. J. Buschow, 
J. Appl. Phys. 78, 5568 (1995). 

26 D. Hautot, G. J. Long, F. Grandjean, and O. Isnard, Phys. 
Rev. B 62, 11731 (2000). 

27 Y. Janssen, M. Angst, K. W. Dennis, R. W. McCallum, 
and P. C. Canfield, J. Cryst. Growth 285, 670 (2005). 

28 O. Isnard, S. Miraglia, D. Fruchart, C. Giorgetti, 
S. Pizzini, E. Dartyge, G. Krill, and J. P. Kappler, Phys. 
Rev. B 49, 15692 (1994). 

29 O. Isnard, S. Miraglia, D. Fruchart, C. Giorgetti, E. Dar- 
tyge, and G. Krill, J. Phys. Condens. Matter 8, 2437 



18 



(1996) . 

30 O. Isnard, S. Miraglia, C. Giorgetti, E. Dartyge, G. Krill, 
and D. Fruchart, J. Alloys Compd. 262-263, 198 (1997). 

31 D. Vandormael, F. Grandjean, V. Briois, D. P. Middleton, 
K. H. J. Buschow, and G. J. Long, Phys. Rev. B 56, 6100 

(1997) . 

32 A. Kreyssig, S. Chang, Y. Janssen, J. .-W. Kim, S. Nandi, 
J. Q. Yan, L. Tan, R. J. McQueeney, P. C. Canfield, and 
A. Goldman, Submitted to Phys. Rev. B (2006). 

33 P. C. Canfield and Z. Fisk, Philos. Mag. B 65, 1117 (1992). 

34 P. C. Canfield and I. R. Fisher, J. Cryst. Growth 225, 155 
(2001). 

35 K. Koyama, H. Fujii, and P. C. Canfield, Physica B 226, 
363 (1996). 

36 L. H. Lewis, J. Y. Wang, and P. C. Canfield, J. Appl. Phys 
83, 6843 (1998). 

37 Materials Preparation Center, Ames Laboratory US-DOE, 
Ames IA 50011, http://www.mpc.ameslab.gov 

38 H. Okamoto, Desk Handbook Phase diagrams for Binary 
Alloys (ASM International, Materials Park, OH, 2000). 

39 T. C. Ozawa and S. J. Kang, J. Appl. Cryst. 37, 679 
(2004). 

40 G. Calestani, N. Magnani, A. Paoluzi, L. Pareti, and 

C. Rizzoli, Phys. Rev. B 68, 054424 (2003). 

41 O. Isnard, S. Miraglia, J. L. Soubeyroux, D. Fruchart, and 
A. Stergiou, J. Less-Common Met. 162, 273 (1990). 

42 R. D. Shannon, Acta Cryst. A32, 751 (1976). 

43 E. Girt and Z. Altounian, Phys. Rev. B 57, 5711 (1998). 

44 B. Saje, S. Kobe-Besenicar, Z. Samardzija, A. E. Platts, 

D. Kolar, and I. R. Harris, J. Magn. Magn. Mater. 146, 
L251 (1995). 

45 D. H. Dennison, M. J. Tschetter, and K. A. Gschneidner, 



Jr., J. Less-Common Met. 10, 108 (1966). 

46 M. Ellerby, K. A. McEwen, E. Bauer, R. Hauser, and 
J. Jensen, Phys. Rev. B 61, 6790 (2000). 

47 M. Ellerby, K. A. McEwen, and J. Jensen, Phys. Rev. B 
57, 8416 (1998). 

48 J. M. Ziman, Electrons and Phonons, The Theory of Trans- 
port Phenomena in Solids (Oxford University Press, Lon- 
don, 1963). 

49 H. Fukuda, Y. Janssen, H. Fujii, T. Ekino, and Y. Morii, 
Nippon Oyo Jiki Gakkaishi 23, 108 (1999). 

50 R. A. Ribeiro, S. L. Bud'ko, and P. C. Canfield, J. Magn. 
Magn. Mater. 267, 216 (2003). 

51 E. Morosan, S. L. Budko, P. C. Canfield, M. S. Torikachvili, 
and A. H. Lacerda, J. Magn. Magn. Mater. 277, 298 
(2004). 

52 M. Angst, A. Kreyssig, Y. Janssen, J.-W. Kim, L. Tan, 
D. Wermeille, Y. Mozharivskyj, A. Kracher, A. Goldman, 
and P. Canfield, Phys. Rev. B 72, 17440 (2005). 

53 M. E. Fisher, Phil. Mag. 7, 1731 (1962). 

54 W. P. Wolf and A. F. G. Wyatt, Phys. Rev. Lett. 13, 368 
(1964). 

55 L. J. de Jongh and A. R. Miedema, Adv. Phys. 23, 1 
(1974). 

56 E. Stryjewski and N. Giordano, Adv. Phys. 26, 487 (1977). 

57 T. Moriya and K. Usami, Solid State Commun. 23, 935 
(1977). 

58 L. D. Landau, Zh. Eksp. Teor. Fiz. (Sov. Phys. JETP) 7, 
19 (1937). 

59 S. K. Yip, T. Li, and P. Kumar, Phys. Rev. B 43, 2742 
(1991). 

60 H. G. M. Duijn, PhD. Thesis (Universiteit van Amster- 
dam, Amsterdam, 2000). 



